Tuning a Josephson junction through a quantum critical point 
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We tune the barrier of a Josephson junction through a zero-temperature metal-insulator transition 
and study the thermodynamic behavior of the junction in the proximity of the quantum-critical 
point. We examine a short-coherence-length superconductor and a barrier (that is described by a 
Falicov-Kimball model) using the local approximation and dynamical mean-field theory. The inho- 
mogeneous system is self-consistently solved by performing a Fourier transformation in the planar 
momentum and exactly inverting the remaining one-dimensional matrix with the renormalized per- 
turbation expansion. Our results show a delicate interplay between oscillations on the scale of the 
Fermi wavelength and pair-field correlations on the scale of the coherence length, variations in the 
current-phase relationship, and dramatic changes in the characteristic voltage as a function of the 
barrier thickness or correlation strength (which can lead to an "intrinsic" pinhole effect). 
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I. INTRODUCTION 

The theoretical understanding of Josephson junctians 
has progressed dramatically over the past four decadesd'H. 
Recent advancesci have been fostered hv, nanofabrication 
of superconducting mesoscopic devicesJJ which, together 
with high-temperature superconductor junctions, have 
revived interest in the transport properties of supercon- 
ductors weakly coupled through a normal region. The in- 
terplay between phase-coherent electron propagation in 
the normal region and macroscopic phase coherence of 
Cooper pairs in superconductors generates novel quan- 
tum interference phenomena since the proximity effect in 
such systepjis is mediated by a phase-coherent Andreev 
reflectionB'Q. However, little attention has been paid to 
quantum effects on transport arising from many-body 
correlations in the barrier separating the superconduc- 
tors. Such junctions are frequently encountered in high- 
Tc systems where both superconducting electrodes and 
the riormal region are highly correlated electronic sys- 
tems .13 

Low-Tc junctions have large superconducting coher- 
ence lengths, and effects on the scale of the Fermi wave- 
length can usually be averaged over to accurately de- 
scribe junctions by a quasiclassical (single-particle) ap- 
proach. As the coherence length of the superconduc- 
tor becomes smaller and smaller (as in high- Tc junc- 
tions) one can no longer ignore the interplay between 
oscillations brought on by the Fermi surface and those 
due to pair-field correlations. In addition, as junction 
sizes are made smaller and smaller, the barrier needs 
to be tuned close to the metal-insulator transition in 
order to maintain a large characteristic voltage (where 
properties of-a Josephson junction have been thought to 
be optimized^). The conventional proximity-effect the- 
ory cannot account for supercurrent transport in junc- 
tions where the barrier approaches a metal-insulator 
transitionQ. Therefore, these junctions must be described 
in a full many-body approach that can properly account 
for the change in character of the quantum mechani- 



cal system as the correlations drive a metal-insulator 
transition. The standard single-particle approaches, like 
theJjull quantum transport theories (scattming formal- 
ismQij'El, and Green's-function techniquesll30)|-er tradi- 
tional quasiclassical Green's-function methodalj are in- 
adequate for this purpose (in general, the quasiclassical 
approaches do not require a quasiparticle assumptionli3, 
but the usual quasiclassical Green's function, employed 
in nonuniform superconductivity problems, can be ex- 
panded in temis of Andreev quasiparticle eigenfunctions 
and energiesH). ■— ■ 
Recent progress in the dynamical mean-field theoryll3 
has shown how to genetaJize the local approximationj-to 
inhomogeneous systemgi3 and to Josephson junctionsEZI. 
Here we utilize this formalism to examine what happens 
as the barrier material is tuned through a quantum- 
critical transition where the single-particle density of 
states is suppressed to zero and a correlated metal- 
insulator transition occurs. We find that in this region of 
phase space, it is important to include self-consistency 
effects and-, many- body effects. The simple analytical 
treatmentsO of Josephson junctions rely on the usage of 
rigid boundary conditions,a i.e., a step function model 
for the pair potential at a normal-superconductor inter- 
face. This is justified in narrow junctions (barrier width 
smaller than the bulk coherence length ^o) where the ef- 
fect of the constriction induced by the narrow barrier 
on the order parameter of the much wider supercon- 
ductors is "geometrically diluted" , or in wide junctions 
with high resistivity barriers (in both cases the critical 
current of the junction is much smaller thaa the bulk 
critical current of the superconducting leadaiS). On the 
other hand, a self-consistent solution for the variation of 
the order parameter A(a::) [i.e., pair-correlation function 
F{x) = —A{x)/U{x) with U{x) the site-dependent inter- 
action strength] induced by the current fiow or geometry, 
not only ensures current conservation and allows pae to 
find the critical current in an arbitrary geometry,ll3 but 
is unavoidable in situations where the proximity effect 
induces appreciable superconductivity in the normal re- 
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gion^pOr when the thickness of the weak Hnk is com- 
parablellj tp-/o- Thus, the microscopic self-consistent 
calculationgl3 reveal a variation of A on length scales 
(like Xp, the Fermi wavelength) smaller than (which 
is also of importance in high- Tc junctions where the qua- 
siclassical approximation,E3 i^o ^ ^f, docs not hold). 
Our junctions are wide, and even in the tunneling limit 
(i.e., with a correlated insulator barrier), they require a 
self-consistent treatment because the many-body effects 
prevent a description in terms of simple phenomenologi- 
cal parameters (like the barrier transparency) . Both self- 
consistency effects and many-body correlations are auto- 
matically included via the dynamical mean-field theory. 

Our results should shed light on high- Tc superconduc- 
tors even though we are restricting ourselves to s-wave 
symmetry order parameters. This is because the high- 
Tc superconductors have short coherence lengths (on the 
order of a few lattice spacings) and have barrier materi- 
als [either from grain boundaries, ion-damage, or doping 
(such as Co-doping)] that are correlated and lie close to 
the Mott metal-insulator transition!]. Our examination 
of s-wave superconductors in this limit illuminates this 
new physical regime without adding the complicated ge- 
ometrical effects that arise from d-wave order parameters 
(which will be investigated in a future study). 

In Section II we briefly describe the formalism that 
is used in our computational techniques. Section III 
contains results on tuning through the quantum criti- 
cal point by increasing the correlation energy at a fixed 
barrier thickness. We examine four cases: (i) thin bar- 
rier; (ii) bilayer barrier; (iii) barrier on the order of the 
bulk coherence length; and (iv) thick barrier. In Section 
IV we tune the metal-insulator transition by increasing 
the thickness of the barrier at fixed correlation energy. 
We examine a weakly-correlated metal barrier, a strongly 
correlated metal, and a Mott insulator, finding deviations 
from quasiclassical results for the correlated insulator. 
Our conclusions are presented in Section V. 



II. FORMALISM 

The computations require a self-consistent calculation 
of the properties of a Josephson junction within a many- 
body formalism. To start, we need to have a solution of 
the bulk superconductor, which will provide the "bound- 
ary condition" for the simulations in the bulk bound- 
aries of the junction. The bulk problem can be solved 
directly in both the absence of a supercurrent and in 
the presence of a supercurrent (where there is a uni- 
form variation in the phase of the superconducting or- 
der parameter). The uniform bulk solution is then em- 
ployed to provide the boundary conditions for the junc- 
tion beyond the region where we determine properties 
self-consistently. The inhomogeneous problem consists 
of N self-consistent planes embedded in the bulk super- 
conductor on each side (see Fig. |l|). The self consistent 



region consists of a sandwich of A'';, barrier planes sur- 
rounded by Nsc planes on each side N = Nb + 2Nsc (the 
word "barrier" is used since the material through which 
the weak link between superconductors is made will have 
its properties tuned from a metal to an insulator) . In our 
solutions we choose Nsc — 30 and TVf, ranges from 1 to 
80. Since the coherence length of the superconductor is 
^0 = hvp/iT^^) ~ 4a (with a the lattice spacing) the self- 
consistent superconducting region is approximately eight 
times the bulk coherence length, which we find to be suf- 
ficient for our calculations. This approach is useful, since 
it does not require us to make any assumptions about the 
boundary conditions at the interface between the bar- 
rier and the superconductor, since they are determined 
self-consistently. The approximation is the presence of 
a (typically small) discontinuity in the supercurrent at 
the bulk superconductor-self-consistent-superconductor 
interface. We have found that the superconducting or- 
der has always healed to its bulk value at that point, 
but sometimes there can be a jump of the supercon- 
ducting phase when one nears the critical current of the 
junction. This discontinuity in the phase (correspond- 
ing to a breakdown of current conservation at the bulk 
superconductor-simulated superconductor interface) can 
become large for thick insulating barriers or when one 
lies on the decreasing current side of the current-phase 
diagram (see below). 
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FIG. 1. Microscopic stacked planar geometry of a Joseph- 
son junction. The sandwich of iV = 10 planes; Nsc = 4 
superconducting planes coupled to a bulk superconductor on 
the left and Nt = 2 barrier planes on the right, followed by 
a further Nsc = 4 superconducting planes coupled to another 
bulk superconductor on the right. The junction is allowed 
to have spatial inhomogeneity only within the A'^ modeled 
planes, but the calculations are for an infinite system. In our 
calculations we always take Nsc = 30 and Nb ranges from 1 
to 80. 

We simulate an inhomogeneous sys- 
tem of stacked square lattice planes that correspond to 
the superconductor-barrier-superconductor sandwich of a 
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Josephson junction. A lattice site corresponds to a unit 
cell (which we normally picture as being one atomic site) , 
and we assume a tight-binding picture with the same hop- 
ping integral Uj between atomic sites within a plane and 
atomic sites between the planes. This description implies 
that we are assuming the "bare" kinetic energy of the su- 
perconductor and the barrier are identical (note that the 
renormalized density of states can be very different, es- 
pecially when the barrier is a correlated insulator). Such 
a condition is not necessary in this formalism, but we 
include it for simplicity, since it reduces the number of 
parameters that are varied in the junction. The supefc. 
conductor is described by an attractive Hubbard modelEll 
in the Hartree-Fock approximation. This is equivalent to 
the conventional Bardeen-Cooper-Schripffeir3 (or more 
accurately, the Boguliubov-de GennesE^ [which involves 
a nonconstant density of states due to the tight-binding 
approach]) description, except in this case the energy 
cutoff is determined by the electronic bandwidth rather 
than the phonon frequency. In fact, the attractive Hub- 
bard model offers richer behavior (that is not employed in 
this contribution) showing a crossover to preformed pairs. 
By including higher-order processes in U, tiirough a T- 
matrixEJ or a dynamical niean-field theoryCJ approach, 
one can study the crossoveiEa from BCS superconductiv- 
ity, where pair formation and condensation occur at 
to preformed pairs that condense at a lower temperature 
(this should be important in short-coherence-length su- 
perconductors like the high- Tc materials]. The barrier 
is described by a Falicov-Kimball modelED at half filling. 
This model has two kinds of particles: (i) mobile elec- 
trons and (ii) static ions. The average concentration of 
electrons is one per site and the average concentration of 
ions is one-half per site. When an electron and an ion oc- 
cupy the same lattice site, there is a Coulomb attraction 
Ufk between them. One can view this system as a binary 
alloy of A and B ions at 50% concentration with Ufk 
being the difference in site energy between the A and B 
ionic sites (the off-diagonal energy is assumed to be the 
same for the A and B ions). The many-body problem is 
solved by taking an annealed average and is essentially 
the simplest disorder problem (and the simplest many- 
body problem) . It undergoes a metal- insulator transition 
in the bulkEa (see below) which is why we adopt it for 
study here. 

The Hamiltonian for the Josephson junction is then 

+Y^urcic.Jw.-l), (1) 

where cf^. {cia ) creates (destroys) an electron of spin a at 
site i, tij is the hopping integral between nearest neigh- 
bor sites i and j (we measure energies in units of i), 
Ui = —2 is the attractive Hubbard interaction for sites 
within the superconducting planes, Uf'^ is the Falicov- 



Kimball interaction for planes within the barrier, and Wi 
is a classical variable that equals 1 if an A ion occupies 
site i and is zero if a i? ion occupies site i. A chemical 
potential fi is employed to determine the filling. Since we 
work at half filling for both the superconductor and the 
barrier, we have ^ — 0. Note that if Ui ~ Uf^ — for 
all lattice sites, the Hamiltonian describes tight-binding 
electrons on a simple cubic lattice. 

The superconducting regions are described by an at- 
tractive Hubbard model with Ui = —2 and Wi = for 
all superconductor sites. The homogeneous bulk super- 
conductor has a transition temperature Tc — 0.11 and 
a zero-temperature order parameter A — 0.198. This 
yields a standard BCS gap ratio 2A/ {ksTc) ~ 3.6 and a 
coherence length = fivp/{TrA) that ranges from 3.5a 
to 4.3a depending on whether we average the absolute 
value of Wi? over the Fermi surface or take the root-mean- 
square oivp (a cubic lattice at half-filling has a direction- 
dependent Fermi velocity); a fit of the decay of the super- 
conducting order^s it is disturbed at the superconductor- 
barrier interfacetZl gives ~ 3.7a. The bulk critical cur- 
rent per unit area is /c,buik = 0.0289(2ei)/(?ia^). The 
value of our bulk critical current density is slightly higher 
than the one determined by a Landau depairing velocity 
Vd — A/hkp (jc.buik — envd, where the density of parti- 
cles is n = fc|,/27r^, assuming a spherical Fermi surface) 
because of the possibility to have gapless superconduc- 
tivity in tiiKce dimensions at superfluid velocities slightly 
exceedingcj Vd (note that kp is direction-dependent for a 
cubic lattice at half-filling). Calculations on our junction 
are performed at a temperature of T = 0.01, which is ef- 
fectively at the zero-temperature limit {T/Tc « 0.09) for 
the superconducting properties. The barrier region is de- 
scribed by a half-filled Falicov-Kimball model in the sym- 
metric limit of {wi) = 0.5. In the bulk, this barrier un- 
dergoes a metal-insulator transition at UpK ~ 4.9 (since 
the bandwidth of the simple cubic lattice is 12, the metal- 
insulator transition occurs when U fk is on the order of 
one- half of the bandwidth) . This is illustrated in Fig. ||(a) 
where we show the single-particle density of states for a 
bulk barrier as a function of UpK ■ The density-of states 
for this model is independent of temperatureEa. Since 
the system is not a Fermi liquid for nonzero Upx, one 
can see the density of states first develops a pseudogap 
and then is suppressed entirely to zero as the correlations 
are increased and it becomes a correlated insulator. The 
opening of the gap is continuous. In Fig. ||(b) we show 
the imaginary part of the local self energy at low-iCner- 
gies. This result is also temperature independentc3. We 
see that the curvature of the self energy has the wrong 
sign in the metallic regime (which is one reason why it is 
not a Fermi liquid) and that it diverges (and becomes a 
delta function) as the system crosses over into the insulat- 
ing phase. This occurs because the self energy develops 
a pole at zero energy in the insulating phase. Such be- 
havior can only be seen in a many-body treatment of the 
system. 

The inhomogeneous system is solved by employing the 
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FIG. 2. (a) Electronic density of states for the bulk barrier 
(simple-cubic lattice) described by the Falicov-Kimball model 
in the local approximation. The value of Ufk ranges from 1 
to 7 in steps of 1. As Ufk increases the density of states 
first develops a pseudogap and then a real gap. (b) Absolute 
value of the imaginary part of the local retarded self energy 
for low frequency on the real axis. See how the curvature 
has the wrong sign for a Fermi liquid and how the imagi- 
nary part diverges at zero frequency as one goes through the 
quantum-critical point and a pole develops in the self energy. 

matrix formulation of NambuS for the Green's function 
G{ri,rj,iuJn) for two lattice sites and Vj at the Mat- 
subara frequency iun = i7rT(2n -|- 1), 

^/ ^ ■ , ^ f G{r^,rj,iujn) F{r,,rj,iujn) \ 

(2) 

and the corresponding local self energy, 

E(r„zc.„)=('5,'"-*^^") t^:;'*'^"\V (3) 

The diagonal and off-diagonal Green's functions are de- 
fined respectively as: 

G{r„rj,iujn) = - dTexp(zcj„T)(T^CjCT(r)cj^(0)), (4) 
Jo 
rP 

F(rj,rj,zw„) = - / dTexp(zw„T) (T^Cji(t)cji(O)) , (5) 
Jo 

where T,- denotes time-ordering in r and (3 = l/T. 

The self energies and Green's functions are coupled 
together through Dyson's equation, 

G{ri,rj,iuJn) = G'°(ri, r^, icj„) 

+ ^G°(ri,r/,iw„)S(r;,iw„)G(r/,rj,iw„), (6) 



where we have included the local approximation for the 
self energy, S(ri,rj,jw„) = S(ri, ia;„)(5ij . The non- 
interacting Green's function, G°(ri, , iw^) is diagonal 
in Nambu space, with upper diagonal component given 
by: 

G" (r, , r, , zc^„) = / d^k-^^^^^-^ . (7) 

J IWn -I- /i - £k 

We emphasize that G" is the non-interacting Green's 
function and is not the effective medium of an equiva- 
lent atomic problem. 

Details of the |-Gpmputational scheme have been de- 
scribed elsewheretll. Here we simply summarize the algo- 
rithm. The junction is inhomogeneous in the z-direction 
only, since it has translational symmetry within each 
plane. The algorithm begins by converting the three- 
dimensional system to a quasi-one-dimeusional system 
using the method of Potthoff and NoltingE£l. We perform 
a Fourier transformation within each plane to determine 
the mixed-basis Green's function [defined in terms of two- 
dimensional momenta {kx and ky) and the z-coordinate 
of the plane] under the assumption that the electronic self 
energy is local (but can vary from plane to plane). For 
each momentum in the two-dimensional Brillouin zone, 
we have a one-dimensional problem with a sparse ma- 
trix, since the only coupling between planes is due to the 
hopping to each neighboring plane. The infinite "tridiag- 
onal" matrix can be invertedJpy employing the renormal- 
ized perturbation expansionEil, which calculates both the 
single plane and the nearest neighbor Green's functions. 
A final summation over the two-dimensional momenta 
produces the local Green's function and the Green's func- 
tion for propagation from one plane to its neighboring 
plane. The dynamical mean-field theory is then employed 
to calculate the local self energy from the local Green's 
function and then the local Green's function is calculated 
from inverting the quasi-one-dimensional matrix. These 
two steps are repeated until the Green's functions have 
converged to a fixed point. At the fixed point, we have 
a self-consistent solution of the inhomogeneous problem 
that allows for nonuniform variations in both the pair- 
field correlations (or equivalently the superconducting or- 
der parameter) and in the phase. One important consis- 
tency check is total current conservation at each plane 
in the self-consistent region. All calculations conserve 
current except in extreme cases for thick insulating bar- 
riers (see below). But there can be discontinuities in the 
current at the bulk-superconductor-self-consistent super- 
conductor interface (since this is far from the Josephson 
junction, it has a negligible effect on the results). This 
computational algorithm is a generalization of the con- 
ventional Boguliubov-de Gennes approach to allow for 
correlations within the barrier. 

This algorithm can be performed for the normal state 
or for the superconducting state and can be performed 
on the imaginary or real frequency axes. We work on the 
real axis in order to calculate the normal state resistance. 
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Since we have a many-body system, we must use Kubo's 
formula for the con|ductivity. Details for this calculation 
appeared elsewheretZl. Our formalism calculates the con- 
ductivity by neglecting vertex corrections and evaluating 
the simple bubble diagram (which becomes exact in the 
infinite-dimensional limitcj) . 



III. TUNING THE CORRELATION STRENGTH 
THROUGH A METAL-INSULATOR 
TRANSITION 

We begin by presenting results for a fixed barrier thick- 
ness, and vary the Falicov-Kimball coupling strength. 
We study four different systems: (i) a thin barrier with 
iVb = 1; (ii) a bilayer barrier with Nj, = 2; (iii) a bar- 
rier on the order of the bulk superconducting coherence 
length Nf, — 5; and (iv) a thick barrier A^;, = 20. 



A. Thin barrier [Nb = 1) 

A single plane barrier must be in the very strong in- 
sulating limit before it can severely affect the transport 
perpendicular to the plane. Hence we expect to have 
to increase the Falicov-Kimball interaction to be much 
larger than 5 before the junction starts to display "insu- 
lating" behavior. Similarly, we expect the critical current 
to be close to the bulk critical current, because the plane 
is so thin (at least for metallic barriers). In this regime, 
self-consistency— is critical in determining the properties 
of the junctior£3. 

We begin by examining the proximity effect within the 
junction. Since the Hubbard attraction is zero within 
the barrier, the superconducting gap A, which is pro- 
portional to the Hubbard attraction, identically vanishes 
there. But we can still examine the superconducting pair- 
field correlations by plotting the anomalous average at 
equal times F{t = 0+). This Green's function is contin- 
uous as one passes through the superconductor-barrier 
interface. We show i^(0+) in Fig. ^ Notice how the cor- 
related metal {Ufk < 2.5) appears just as we expect it 
to: the superconductivity is smoothly depressed as we ap- 
proach the barrier and then decreases within the barrier 
as correlations increase. As the bulk barrier enters the 
pseudogap regime (2.5 < Ufk < 5) we see small oscil- 
lations appear in the superconductor, and the supercon- 
ductivity continues to be depressed within the barrier. 
In the correlated insulator regime (5 < Ufk < 8) the 
oscillations continue to grow and the superconductivity 
within the barrier is small, but rather insensitive to Ufk- 
In the strong insulating regime {Ufk > 8), we find that 
the oscillations become large and the superconductivity 
eventually becomes enhanced within the barrieri We be- 
lieve that the oscillations and this enhancement of the 
anomalous average are arising from a surface effect of 



the superconducting half-planes — each half-plane devel- 
ops oscillations near the surface (as the barrier becomes 
more insulating). 
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FIG. 3. Proximity effect for a thin barrier {Nb — 1). The 
anomalous average is plotted versus plane number (the in- 
sulating barrier lies at plane 31). The numbers indicate the 
value of Ufk (1,2,3,4,6,8,12,16, and 20); the anomalous aver- 
age monotonically increases with Ufk in the range between 
planes 29 and 30. Note how oscillations develop as the correla- 
tions increase until the superconductivity ultimately becomes 
enhanced in the barrier for the strong insulator. 

In Fig. ^, we show plots of (a) the current-phase rela- 
tion as well as (b) the normalized relation I{9)/Ic. The 
phase difference across the junction is defined as the to- 
tal phase across the barrier plane. Since there is only 
one plane, and the majority of the phase jump occurs 
at the barrier plane, we find that the Nf, = 1 Joseph- 
son junction has significant phase change over the super- 
conducting region, since we must define the barrier to 
begin at a distance halfway between the last supercon- 
ducting plane (on the left) and the barrier plane and end 
halfway between the barrier plane and the first super- 
conducting plane (on the right). This result arises from 
the fact that the barrier plane is so thin and because we 
have discretized real space to correspond to the atomic 
unit cells. Note how the shape of the current-phase re- 
lation is far from sinusoidal when the critical current is 
close to the bulk critical current. As the barrier becomes 
more insulating, the critical current decreases and ap- 
proaches a more sinusoidal shape. Because of the lattice 
nature of the model, the current-phase relation actually 
approaches I{9) = U sm{26) for a single-plane junction 
(note that this is an artifact of the coarse-graining in 
real space that maps the junction onto a lattice). More 
metallic barriers have the maximum of the current occur 
at a phase difference much smaller thaii.7r/4 as predicted 
to occur in seli-consistent calculationsEa, and seen in our 
previous workll3. 
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FIG. 4. Current-phase relation for a single-plane bar- 
rier {Nb — 1). (a) I{9) plotted for various values of Ufk 
(1,2,3,4,6,8,12,16, and 20); the numbers denote the value of 
Ufk- The bulk critical current is shown for comparison as 
the horizontal dashed line. Note how in the weakly corre- 
lated limit the curve is nonsinusoidal, but becomes sinusoidal 
as the correlations increase, (b) Normalized 1(6) /Ic- Note 
how this approaches sin 20 for strong correlation, as discussed 
in the text. [The same values of Ufk are plotted in (b), but 
we don't label the curves here.] 



Finally we show a plot of the Josephson critical cur- 
rent, the normal state resistance i?„, and the charac- 
teristic voltage IcRn for the single plane barrier. The 
plots are on a semilogarithmic scale. The critical cur- 
rent drops by about two orders of magnitude as U fk 
increases from to 20. This occurs even though the 
anomalous average increases in the strongly insulating 
limit! The plot of i?„ shows the expected increase as 
the correlations increase. But even at a large value of 
Ufk — 20, the resistance only increases by about two or- 
ders of magnitude over the noninteracting limit. Finally, 
we show the characteristic voltage IcRn in Fig. |^(c). Its 
value does not change much, but shows a mild optimiza- 
tion for the moderately correlated metallic phase. In the 
metallic limit [/ ^ 0, the /ci?n product approaches the 
product of the bulk critical current times the Sharvin 
resistance, which is 0.287t/e — 1.45A/e. This resuljt is 
different from the clean Kulik-Omelyanchuk limitcj of 
TT A / e because we are treating a different geometry from 
a point contact (which can be described as a "plane" con- 
tact). As the correlations increase within the metallic 
phase, the characteristic voltage peaks for U fk ~ 2 at a 
value somewhat smallp:.than the dirty limit of the Kulik- 
Omelyanchuk formulao for a superconducting point con- 
tact IcRn — 0.667rA/e at T = 0. In fact, there are two 
possibilities for the IcRn product (i.e., critical current) 
of a short contact with diffusive scattering. Namely, in 
the single-particle picture scattering properties of a nor- 



mal region can be described by the universal distribu- 
tion of transparencies D (defined as the distribution of 
eigenvalues of the matrix tt^, where t is the transmis- 
sion matrix connecting incoming to outgoing transverse 
propagajtiiig modesH) given by either the Dorokhov ex- 
pressionN PYyo{D) = (G/2Gq) [D^JT^]-^ (valid for 
most bulk conductors) , or the Schep-Bauer distributionEj 
Psb{D) = {GJ^Gq) [1)3/2 VT^n?]-! (-yj^ijj f-Qj. sub-nm- 
thick barriersL3). Here, Gq — jh is the conductance 
quantum and G = dD P{D)D is the disorder-averaged 
conductance. The total current is found by integrat- 
ing the current carried by a single channel I{D) (with a 
transparency of D) over the distribution function. P{D) 
as shown in the multiple Andreev reflection theorytl. This 
integral, I — dD P{D)I{D), then leads to the follow- 
ing characteristic voltages: IcRn = 0.667rA/e for Pdo(-D) 
and IcRn = 0.6l7rA/e for Psb{D). Wc find that in the 
case of a single-plane barrier made of an FK correlated 
metal, the largest IcRn (obtained for U fk = 2) is slightly 
below the value determined by Psb{D)- This can be at- 
tributed to effects of self-consistency (which always lower 
the critical current because of the depression of the order 
parameter in the superconducting leads due to the prox- 
imity effect), or to the fact that such an interface cannot 
be described by the Schep-Bauer distribution Psb (the 
rationale behind the comparison of our barriers, domi- 
nated by many-body correlations, with conventional re- 
sults relying on the single-particle picture of transport, 
is elaborated further in Sec. 
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As the junction barrier becomes more insulating, the 
characteristic voltage becomes essentially coiislant as ex- 
pected from the Ambegaokar-Baratoff limitO tt A / (2e) 
(dashed line). But the magnitude of the characteris- 
tic voltage is approximately 15% smaller than that pre- 
dicted by them [7rA/(2e) versus our calculated result of 
1.31A/e]. Once again, this small reduction arises from 
the self-consistency for a short-coherence-length super- 
conductor which reduces the gap as one approaches the 
barrier and from the Fermi-surface averaging of the trans- 
port, since the Fermi surface is far from spherical at half- 
filling. 



B. Bilayer (Nt = 2) 

We see similar behavior in the bilayer junctions with 
Njj ~ 2. The correlation strength needed to make an 
insulating barrier is smaller here, because the barrier is 
thicker. In Fig. ^ we plot the anomalous average as a 
function of the plane number. The planes numbered 31 
and 32 are where the barrier lies. The behavior is like 
that seen in the thin barrier case — the correlated metal 
{Ufk < 2.5) and pseudogap regions (2.5 < Ufk < 5) 
are similar. The correlated insulator regime, where the 
oscillations in the anomalous average increase, but its 
value within the barrier is rather insensitive to Ufk (5 < 
Ufk < 7) and the strong insulating regime {Ufk > 7) 
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shows an enhancement of the anomalous average within 
the barrier at an even larger value of U fk (starting at 
Ufk ~ 14). 



0.16 



10 



£^ 10-^ 



o 10" 

> 10 



-6 



:^ 10" 

^ 10"^ 

10^ 

10° 



1 1 1 1 1 1 1 


1 1 1 1 1 


(a)- 







. . 1 , , . . 1 


5 


10 


15 20 


1 1 1 1 1 1 1 

:(b) 


1 1 1 1 1 


1 1 1 1 







. . 1 , , . . 1 


5 


10 


15 20 


_ 1 1 1 1 1 1 1 

(c) 


1 1 1 1 1 


1 1 1 1 1 1 1 1 _ 


^ • • • • 


1 , , . . 1 



5 10 15 

Interaction strength Up|^ 



20 



Freericks et al, Phys. Rev. B, Figure 5 
FIG. 5. (a) Critical current, (b) normal state resistance, 
and (c) characteristic voltage of the Josephson junctions as 
a function of the Falicov-Kimball interaction within the bar- 
rier. The circular symbols are for A^;, = 1 and the squares 
are for Ni, = 2. Note how the critical current decreases and 
the junction resistance increases as expected, and how the 
characteristic voltage does not depend too strongly on the 
correlation strength. The dependence on correlation strength 
for the bilayer junction is much stronger than for the thin 
junction. The dashed line in (c) is the Ambegaokar-Baratoff 
prediction. 
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Freericks et ol, Phys. Rev. B, Figure 6 

FIG. 6. Anomalous average for differing correlation 
strengths and Nt — 2. Note how similar these results are 
to the thin barrier case — how the anomalous average initially 
drops within the barrier, then oscillations develop, followed 
by an enhancement for the strongest correlation strengths. 

In Fig. 0, we show (a) the current versus phase and 
(b) the renormalized current-phase relation. The current 
is reduced by a factor of about 25% for an equivalent 
value of the correlation strength, and the phase across 
the junction increases by about a factor of two versus 
the single-plane barrier. As the barrier becomes more 
insulating, we recover the expected result that I{9) = 
Ic sin 6 because now all of the phase difference takes place 
over the barrier (in general, the majority of the phase 
difference occurs over the central plane of the barrier). 
More metallic barriers have the maximum in the I{9) 
curve pushed to values of 6 less than 7r/2 as expected 
due to the self-consistency and the proximity to the bulk 
critical current. 

Self-consistency is an important feature of the calcula- 
tions. One might ask whether the self-consistency mod- 
ifies the superconductivity much when current is being 
carried by the junction versus the zero-current case. We 
find that the self-consistency is modified when the bar- 
rier is thin and it is carrying supercurrent. In particular, 
we find that the anomalous average changes as the phase 
gradient for the bulk system increases. For a weakly cor- 
related metal, the largest change is within the barrier, 
where the anomalous average decreases by about 30% as 
the critical current is approached {Ufk = 2). A smaller 
decrease is seen in the superconducting regions close to 
the barrier. In the strongly-correlated insulator regime, 
the behavior is different — the anomalous average within 
the barrier increases by about 10% as one passes through 
the full I{9) curve {Ufk = 12). The change in the 
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Freericks et al, Phys. Rev. B, Figure 7 

Current phase relations for A^;, — 2. (a) Cur- 
rent-phase relation. Note how the current is typically re- 
duced by about 25% for an equivalent value of Ufk from the 
thin barrier result and how the phase difference is approx- 
imately doubled across the barrier. Values of Ufk include 
1,2,3,4,5,6,8,10,12, and 16 (labels included for some of the 
curves), (b) Normalized current-phase relation. Note how 
the expected sinusoidal dependence enters for insulating bar- 
riers. 



superconducting region is much smaller (less than 1%). 
This modification of the self-consistency when the junc- 
tion carries current becomes less important as the barrier 
thickness increases — in the single-plane barrier case, we 
found the anomalous average within the barrier increases 
by more than a factor of 2 (for Ufk — 12). 

The critical current, normal state resistance and char- 
acteristic voltage appear in Fig. ^. The critical cur- 
rent decreases much more rapidly for the bilayer than 
for the single-layer junction. The normal state resis- 
tance increases more rapidly as well, since the bilayer 
has a resistance that is much more than two times the 
Nb = 1 resistance in the strongly insulating limit. The 
characteristic voltage is quite interesting, because it has 
nonmonotonic behavior. There is a weak maximum for 
the moderately correlated metal (near Ufk = 3), but in 
the insulating region the voltage increases linearly with 
the correlation strength, attaining values more than 40% 
higher than the Ambegaokar-Baratoff limit. This is quite 
different from what we expect — a constant characteristic 
voltage — and the characteristic voltage shows no sign of 
saturating even at a correlation energy of Ufk — 20! 



C. Moderately thick barrier {Nt, = 5) 

The barrier region (with Ni, = 5) is chosen to be 
slightly thicker than the bulk coherence length ~ 4a. 
Here we examine properties of a junction outside of the 
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Freericks et al, Phys. Rev. B, Figure 8 

FIG. 8. Anomalous average plotted versus plane num- 
ber for the Nb = 5 junction. The values of Ufk chosen are 
1,2,3,4,5,6,7,8, and 10. Note how the large oscillations are 
now separated from each other and are clearly tied to the 
sup erconduct or- barrier int erf ace. 



regime of most analytic approximations. In Fig. we 
plot the anomalous average versus plane number (the 
barrier lies at planes numbered 31 to 35). These results 
are similar to those seen before. In the weak correlation 
regime {Ufk < 2.5), the anomalous average is a smooth 
function that decreases as the correlations increase. Os- 
cillations begin to develop for 2.5 < Ufk < 4.5, but 
the anomalous average continues to decrease within the 
barrier. As the correlations increase further, Ufk > 5, 
the anomalous average first increases at the center of the 
barrier, then a two-peak structure emerges, which has a 
large amplitude oscillation and a minimum at the cen- 
tral plane of the barrier. We can see clearly here that 
the oscillatory behavior seen in the previous cases is aris- 
ing from effects occurring at the superconductor-barrier 
interface as the barrier is tuned through the quantum- 
critical point (this is further confirmed with the Ni, = 20 
data below). 

The current-phase relation is similar to those seen pre- 
viously, and will not be shown here. The normalized 
current-phase relation, is plotted in Fig. |9| This re- 
sult is quite interesting. In the weakly correlated regime 
Ufk < 2.5, the maximum of the current-phase relation 
occurs at a phase smaller than tt/2 as expected for a thin 
metallic barrier. As the correlations increase, the maxi- 
mum first overshoots 7r/2 {Ufk = 3, 4), and then returns 
to its expected location at 7r/2 for Ufk > 5. There is a 
delicate interplay between the strength of the correlations 
and the location of the maximum of the current-phase re- 
lation. 

The critical current, normal-state resistance, and char- 
acteristic voltage are plotted in Fig. |l^ for Nb — 5 (dia- 
mond). It is difficult to locate the metal- insulator tran- 
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Freerlcks et al, Phys. Rev. B, Figure 9 

FIG. 9. Normalized current-phase relation. Note how the 
maximum lies at an angle less than 7r/2 for Ufk = 1,2, in- 
creases to a value larger than 7r/2 for moderate correlations 
Ufk =3,4 and then settles down to ■w/2 for larger correlar- 
tions U FK > 5 (curves with a dotted line) . 

sition from the iV(, = 5 critical current data (except by 
focusing on the inflection point), but the transition is 
clear in the resistance, which has a sharp increase in the 
range from Ufk = 5 to 6. The characteristic voltage 
has striking behavior. Starting at a value about 20% 
less than the Ambcgaokar-BaratofF limit in the metallic 
regime, the voltage initially decreases with correlation 
strength, then has a sharp increase (by over 100%) at the 
metal-insulator transition, reaching a maximum almost 
40% higher than the Ambegaokar-BaratofF result, until it 
finally starts to decrease as correlations increase further, 
continuing to decrease at the largest value of correlations 
where we performed calculations. Hence, junctions in 
this regime do see an enhancement of the characteris- 
tic voltage on the insulating side of the metal-insulator 
transition. This behavior is quite complex! 

As the barrier size is increased to be on the order of 
the superconducting coherence length we see that there 
is an interplay of a number of diff'erent things: as the 
barrier becomes more insulating, oscillations develop in 
the anomalous average that are pinned to lie near the 
superconductor-barrier interface; the supercurrcnt de- 
pends critically on self-consistency for metallic barriers, 
but as correlations increase, there is an overshoot and 
the maximum of the current-phase relation occurs above 
7r/2 until it settles down to an Ic sin 9 behavior for the 
more insulating barriers; there is little indication of the 
metal-insulator transition in the critical current except 
for the appearance of an inflection point (which occurs 
at the same place that the current-phase relation has 
its maximum move to 7r/2), but the resistance shows a 
clearer picture of the transition. The characteristic volt- 
age is the most interesting. Initially it decreases with 
correlation strength, than has a sharp increase at the 
metal-insulator transition, followed by a maximum and 
a decrease as correlations are increased further. 
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Freericks et al, Phys. Rev. B, Figure 10 

FIG. 10. (a) Critical current, (b) normal state resistance, 
and (c) characteristic voltage of the Ni, — 5 (diamond) junc- 
tion. The metal-insulator transition can be seen in the critical 
current (a) and (more easily) in the resistance (b), as the re- 
gions where the slope of the curves changes most dramatically. 
In the strongly correlated insulating regime, we find the ex- 
ponential decay of the current (and increase of the resistance) 
has a different slope than in the correlated metal regime. The 
characteristic voltage (c) has complex behavior: it first de- 
creases in the metallic regime, then has a sharp increase at the 
metal-insulator transition, followed by a decrease as the corre- 
lations increase further (the Ambegaokar-Baratoff prediction 
is the dashed line). Note how the characteristic voltage is 
maximized just on the insulating side of the metal-insulator 
transition, and how the maximal value is about 40% larger 
than the Ambegaokar-Baratoff prediction. 
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Freericks et al, Phys. Rev. B, Figure 1 1 

FIG. 11. Anomalous average for a thick junction Nt = 20. 
Values of Ufk included are 1,2,3,4,4.5,5,5.5,6,6.5, and 7. Note 
how the metallic regime behaves as expected, but oscillations 
develop in the insulating regime that become huge just inside 
the barrier. 

D. Thick barrier (iVi, = 20) 

The thick barrier junction Nt — 20 behaves in many 
respects like the bulk barrier material. The transition 
from a metal to an insulator occurs at approximately 
Ufk = 5 as in the bulk, and the junction rapidly devel- 
ops oscillations in the anomalous average at UpK ~ 4. 
This is shown in Fig. |ll|. Note how the correlated metal 
regime behaves entirely as expected — the anomalous av- 
erage decreases as one approaches and then enters the 
barrier, but it never gets too small in the metallic regime 
because of the proximity effect. As the correlations in- 
crease, oscillations first develop in the superconductor 
and then move into the interfacial region, penetrating 
about one coherence length into the barrier before they 
are rapidly suppressed within the barrier. We see the 
same phenomenon as in the thinner junctions: in the in- 
sulating regime, the anomalous average can increase to 
above it's bulk value within the barrier, but close to the 
superconductor-insulator interface. As the correlations 
increase, the peak of the anomalous average grows. 

The current-phase relation is essentially sinusoidal in 
these junctions. For the weakly correlated metals, the 
peak of the I{9) curve occurs just above tt/2, but as it 
becomes more insulating, the peak moves downward to- 
ward 7r/2 and the current-phase relation becomes I (9) = 
IcSinO, as expected. There is a computational difficulty 
that enters when we are in the strong-insulator regime for 
a thick-barrier junction. Here the critical current gets ex- 
ponentially small, and the computational algorithm loses 
current conservation through the entire junction (when 
the calculation is halted at a self-consistency error of one 
part in 10^ for the anomalous average at r = 0). Instead, 



we see an exponential decrease of the current from the 
value fixed at the bulk superconductor to the value within 
the barrier. The current is constant within the barrier it- 
self, and the current-phase relation is a nice sine curve, so 
we believe that the critical current found from our algo- 
rithm is accurate, even though, the boundary conditions 
with the bulk are trying to force more current through the 
junction than it can have; i.e., the current discontinuity 
occurs far from the barrier region. This scenario is simi- 
lar to that of Josephson's original analytic schemed, since 
in his case, there is no phase gradient over the supercon- 
ductors, so they carry no current, but there is current in 
the barrier, since there is a phase difference across it. 
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Freericks et al, Phys. Rev. B, Figure 12 

FIG. 12. (a) Critical current, (b) normal state resistance, 
and (c) characteristic voltage of the A^i, = 20 (triangle) junc- 
tions. Note how the thick junction behaves much like the bulk 
material. The metal-insulator transition can be clearly seen 
in the critical current and in the resistance, as the regions 
where the slope of the curves changes most dramatically. In 
the strongly correlated insulating regime, we find the expo- 
nential decay of the current (and increase of the resistance) 
has a different slope than in the correlated metal regime. Note 
that the thick barrier has a sharply suppressed characteristic 
voltage because the decrease in Jc is much sharper than the 
increase in R„. We believe this occurs because the tempera- 
ture dependence of the resistance is strong, even at these low 
temperatures. The dashed line is the Ambegaokar-Baratoff 
prediction. 

The critical current, normal-state resistance, and char- 
acteristic voltage appear in Fig. ^ One can see the 
metal-insulator transition clearly in the Ic curve. Within 
the metal, the critical current has an exponential depen- 
dence on Ufk, within the insulator it has a different 
exponential dependence. In the transition region, it de- 
creases most sharply. We believe the reduction in the rate 
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of decrease in the critical current as Ufk is increased 
into the correlated insulator regime arises in part from 
the fact that the effective junction thickness is thinner 
than the true junction thickness due to the oscillations 
in the anomalous average that develop within the barrier 
region. The resistance shows the expected behavior as 
well. One can clearly see the metal-insulator transition as 
the region where the conductivity changes it's functional 
dependence sharply. Note, however, that the character- 
istic voltage is severely affected by the metal-insulator 
transition, since the decrease in the critical current far 
outweighs the increase in resistance, and the characteris- 
tic voltage decreases by many orders of magnitude as one 
enters the strongly-correlated insulator. One reason why 
this occurs is because the resistance for the correlated in- 
sulator depends strongly on temperature. As T — > 0, the 
resistance becomes very large, but it can be sharply re- 
duced as the temperature increases. Hence, even though 
the critical current is at the zero-temperature limit for 
T/Tc = 0.1, the resistance still has strong temperature 
dependence in this regime, which causes the character- 
istic voltage to be sharply suppressed. Perhaps this be- 
havior plays a role in some junctions that appear to work 
well at low temperatures, but then fail as the tempera- 
ture is increased, because on the insulating side of the 
metal-insulator transition, the resistance has strong tem- 
perature dependence. 

E. Summary of tuning the correlation energy 

We have discovered a number of interesting features 
of Josephson junctions for short-coherence-length super- 
conductors that have their barrier tuned through the 
quantum-critical point of a metal-insulator transition. 
The most striking feature we find is that in the insu- 
lating regime, there are oscillations with a wavelength on 
the order of the Fermi wavelength, that appear at the 
superconductor-barrier interface and decay on the order 
of the coherence length on either side of the interface. 
They can have very large amplitudes (on the order of the 
bulk value of the anomalous average) within the barrier. 
We believe that these oscillations are occurring from a 
"surface" effect intrinsic to the superconductor, and de- 
pending on how close the two interfaces are (determined 
by the thickness of the barrier) they are either indepen- 
dent of each other or can interfere. Note that these re- 
sults differ from those found in-metallic junctions with 
"geometrically diluted" barrierst2l. There, oscillatory be- 
havior was seen even for metallic barriers Ufk = 0, and 
the decay length was much longer, leading to a number of 
cycles before the oscillations are damped. We also found 
interesting results for the current-phase relations. As ex- 
pected, thin junctions typically have Ic occur at a phase 
difference smaller than 7r/2, but in all but the single- 
plane junction, as the correlations increase, the maxi- 
mal Ic occurs at 7r/2 and the curve becomes sinusoidal. 



For thick barriers, we find the maximum occurs larger 
than, but close to 7r/2 for metallic barriers and then mi- 
grates towards 7r/2 as the barrier becomes more insulat- 
ing. Finally, we found new behavior in the characteristic 
voltage of a junction. The characteristic voltage is lim- 
ited in the metallic regime by the bulk critical current 
of the superconductor multiplied by the junction resis- 
tance for a clean barrier (the so-called "planar contact" 
limit). This value is approximately 1.31A/e, which is 
a about 8% smaller than the Ambegaokar-Baratoff re- 
sult for an insulating barrier. As the correlations in- 
crease, Ic decreases to be much below the bulk critical 
current of the junction, and Rn increases. The charac- 
teristic voltage has a rich behavior. For the thin junc- 
tion {Nb = 1) it is maximized in the correlated metallic 
regime, and becomes constant for the insulator. As the 
thickness increases to A*";, = 2, we continue to see a small 
maximum in the metallic regime, but the interesting be- 
havior is that for a wide range of correlation strengths, 
the Ambegaokar-Baratoff result does not hold, and the 
characteristic voltage increases with correlation strength. 
For barrier thicknesses on the order of the correlation 
length (Nh = 5), the behavior is even more complex. 
The voltage initially decreases with correlation strength, 
then has a sharp rise at the metal-insulator transition, 
followed by a maximum for the correlated insulator that 
ultimately decreases as the correlations increase further. 
The Ambegaokar-Baratoff regime doesn't hold here ei- 
ther. Finally in the thick junction regime {Ni, = 20), the 
resistance has a strong dependence on temperature in 
the insulating regime, and even at what is a low temper- 
ature for the superconducting properties, the character- 
istic voltage can decrease significantly as the correlations 
increase. The conclusion that can be drawn from this is 
that one requires a careful tuning of the thickness of the 
barrier, the proximity to the metal-insulator transition, 
and the operating temperature to optimize the proper- 
ties of a junction. This idea is further supported in the 
next section. 



IV. TUNING THE JUNCTION THICKNESS 
THROUGH THE METAL-INSULATOR 
TRANSITION 

Here we present results on tuning the junction from 
the thin to thick barrier at three values of Ufk- (i) 
Ufk = 2 a weakly correlated metal; (ii) Ufk = 4 a 
strongly correlated (pseudogap) metal; and (iii) Ufk — 6 
a correlated insulator. In the correlated metal case, the 
junction can be viewed as an SNS weak link,0 while the 
correlated insulator barrier eventually leads to an SIS 
junction. In the case with Ufk = 2, the normal region is 
a non-Fermi liquid metal that can be described as dirty 
metal (resistivity p„ ~ 240 /ificm with the assumption 
that the lattice constant is 3 A), for Ufk = 4, we get 
a "bad metal" (p„ ~ 2700 fiilcm; such huge resistivj-, 
ties do not necessarily require electronic correlations,E3 
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but are also seen in_model calculations involving-disor- 
dered Feripi liquida^). The early experimentalEB and 
theoreticalcil work on SNS junctions revealed that the 
supercurrent in these structures arises from the prox- 
imity effect: superconducting correlations are generated 
in the normal region with the density of pairs decreas- 
ing exponentially from the SN interface on a scale set 
by normal metal coherence length ^„ = (?iI?/27rfcBT)^/^ 
(here, I? is a classical diffusion constant). The equilib- 
rium current then flows at zero voltage because of the 
overlap of pair-field wave functions from the two super- 
conductors. Recent mesoscopic advances have supple- 
mented this "crudell, ipicture with the analysis of energy- 
resolved quantiticgl3'c3 which become important for phe- 
nomena OH small length scales at low terripesatures and 
voltages.cS The initial theoretical studiescilcd relied on 
Ginzburg-Landau theory (which formally requires T to 
be close to Tc) in the dirty limit (^o ^ ^^ with_£. the 
mean free path) and for long junctions {Ni, ^ Cn)c^- In 
the ensuing approaches, based on quasiclassical Green's 
function formalism, junctions with more general parame- 
ters were described,Bc3 where the proximity effect on the 
superconducting side (i.e., a depression of the order pa^. 
rameter near the SN interface) was taken into accountEJ 
(such effects are treated from the onset in self-consistent 
studies like ours). Thus, the conventional proximity ef- 
fect theories show that the critical current is determined 
primarily by the behavior of the superconducting or- 
der parameter when crossing an SN boundary, while its 
thickness dependence and temperature dependence are 
affected by the way quantum coherence is lost in the nor- 
mal meta,l. However, it is only recently that mesoscopic 
studiesEll have emphasized the importance of the Thou- 
less energy-Exh = KD/N^ = 27rkBT(l/N^ for the prox- 
imity effectll3. Although iSxh is determined by the classi- 
cal diffusion time N"^ /V for a particle to cross the sample, 
it plays a prominent role in various quantum phenom- 
ena encountered in disordered (normal) metal electron 
physics.a In the long junction limit A ^ i?Th, the crit- 
ical current is set by i?Th — according to the, recent qua- 
siclassical (non-self consistent) calculationg£3 eIcRn{T = 
0) 10.82£;Th- In the short junction limit E^xh > A, 
for T ^ 0, the product IcRn is expected to be given by 
the diffusive limit 0.667rA/e of the Kulik-Omelyanchuk 
formulaO, or 1.92A/e in the case of dirty interface with 
Schep-Bauer distribution of transparenciesQ'u (as dis- 
cussed in Sec. HI). The high versus low temperature 
limit is seio by the ratio of ksT and i^xh, or, equiv- 
alently, N\j and since iVf, = ^„ is defined to be the 
length scale at which ksT — E^h- 

While the energy gap A is determined by the (attrac- 
tive) electron-electron interaction in the superconduct- 
ing leads, i?Th is a single-electron concept, and as such 
is not directly applicable to our correlated metal (which 
has no well-defined Landau quasiparticles) . Nevertheless, 
it is a common practice in experimental studiesQ to. ex- 
tract estimates for such "quasiparticle" parameteraD us- 
ing measured values of p„ (even for unusually high resis- 



tance barrier materials, like the underdoped cuprate^), 
and check if the conventional treatment can describe the 
junction. From the formal point of view this corresponds 
to comparing the correlated junction to an SNS weak 
link with a well understood normal metal in the barrier 
region (described in terms of non-interacting quasiparti- 
cles), which has the same resistivity as the real barrier 
material [in fact, most of theoretical treatments oper- 
ate with simple, usually phenomenological, description 
of the barrier region in terms of its transmission prop- 
erties D which reproduce the measured average resis- 
tance 1/Rn = Jq dDP{D)D, as discussed in Sec. III]. 
Therefore, for the sake of comparison of our result with 
standard calculations, we extract a diffusion constant 
V from the Einstein relation l/p„ = 2e'^N{EF)'D [with 
N{Ep) the (single-spin) interacting density of states at 
the Fermi energy] . This is independent of band-structure 
effects (classically V = vplj?), but T) can also be de- 
fined quantum-mechanically ftpm the Kubo formula in 
an exact state representation,^ which then allows one to 
use a diffusivity even when the semiclassical picture of 
the mean free path t breaks down). For the dirty-metal 
case Ufk — 2 we find 'Dup^=2 ~ 2ta'^/h and for the 
bad-metal case we find 2?c/Fif=4 ~ 0.32ta^/?i. The corre- 
sponding normal metal coherence lengths are ^„ « 5.6a 
and ^„ « 2.3a (at T = 0.01) in the former and latter 
case, respectively. 

The critical current is plotted in a semilogarithmic plot 
in Fig. ^ The symbols are the calculated values and the 
dashed lines are a fit to the following form 
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with A a constant, the coherence length in the barrier 
(the symbol £,b is used here to differentiate it from the 
phenomonological ^„ determined from a diffusive metal 
analogy above) , and x an exponent (we only fit the data 
with Nf, > 10, since the thin-plane results are sufficiently 
different from the thick-plane results, and the fits are 
therefore much more accurate; nevertheless, the final 
functional form for the data works well for all barrier 
thicknesses). We find that the fits vary from the ana- 
lytic forms for the thick barrier limit {x = 1). For ex- 
ample we find that the coherence length decreases from 
6 = 6.66 for Ufk = 2, to ^f, = 2.96 for Ufk = 4, 
to = 0.665 for Ufk — 6. Similarly, the exponent 
varies from x — —0.40 for Ufk — 2, to x = —0.45 for 
Ufk = 4, to a; = —0.53 for Ufk = 6. The value for the 
exponent never becomes close to the asymptotic result of 
X = 1 for a thick junction. But the coherence length be- 
haves as expected — as the scattering increases, the coher- 
ence length decreases, becoming very small as the barrier 
goes through the metal-insulator transition and becomes 
a correlated insulator (in fact, our values for are only 
about 20% larger than the estimates for ^„ given above). 
Note that this fitting procedure is not well-defined, since 
we do not have data over many decades of barrier thick- 
nesses and because we can trade-off some effects of the 
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FIG. 13. Critical current versus the barrier thickness for 
(a) a weakly correlated metal {U fk ~ 2, circle) , (b) a strongly 
correlated metal {Ufk = 4, square), and (c) a correlated 
insulator near the quantum-critical point {Ufk ~ 6, triangle). 
Note how the shapes of these curves vary as the correlations 
increase. The dashed curve is a fit of the data for A^t > 5 to 
the form given in Eq. (0). 



fitting by simultaneously changing the exponent and the 
coherence length. But in all cases shovifn, we fit all of the 
data from Ni, = 5 to Ni, = 80 to an accuracy of better 
than 10% for the critical current (the accuracy decreases 
to about 25% for Nb = 1). This fitting scheme with 
nontrivial exponents is definitely more accurate than the 
best fits one can achieve with a; = 1 or a; = 0. We find 
the case with Ufk = 4 to be the toughest case to fit to 
this form (the spread in error is about 10% here, with 
25% error for 7Vb = 1), while the Ufk = 2 is the easiest, 
with a fit for all values A^b = 1 to 80 being accurate to 
5%. 

In Fig. |l^, we plot the characteristic voltage (in units 
of A/e) versus the ratio of the barrier coherence length 
(determined from Fig. ^ and equivalent to the Thou- 
less length) to the barrier thickness. This is our aaa- 
logue of the recent results of a quasiclassical theorycZl, 
which show deviations of the Kulik-Omelyanchuk rela- 
tion for long diffusive junctions. The results for metallic 
junctions Ufk = 2,4, have the same shape as seen in 
the quasiclassical theory, and they nearly scale with each 
other (the scaling could be improved by slightly changing 
the barrier coherence length). The correlated- insulator 
results, Ufk = 6, however, show a different functional 
shape, with the transition from the nearly constant char- 
acteristic voltage to the region where it decreases sharply, 
occurring much more rapidly than in the metallic case 
(and having a small "oscillation" at the "transition"). 



One can be more quantitative in the comparison with 
the quasiclassical predictions: in the kmg-junction limit, 
the characteristic voltage is predictedeH to behave like 



hRn^A'N^' exphiVfc/e,,], 



(9) 



with the coherence length determined from the functional 
dependence of the critical current on N}, in Eq. (^. The 
quasiclassical prediction gives A' « b.AQT)/^^ and x' — 1. 
While we can fit reasonably well to this functional form 
in the regime where £,b/Ni, < 1, we typically find the 
constant A' is about three to five times larger and the 
exponent x' is about a factor of two smaller than the 
quasiclassical predictions for the correlated metal cases 
{Ufk = 2 or 4). The parameters deviate significantly for 
the correlated-insulator phase (where the fitting breaks 
down severely for ^b/^b > 0.1). This shows that the 
correlated-insulator regime cannot be described by the 
conventional quasiclassical approach. There appears to 
be a critical length at which point the characteristic volt- 
age changes from an essentially constant dependence on 
the barrier thickness to a rapidly decreasing dependence 
on the thickness (which is A^;, ~ 7 for Ufk = 6). The 
difference in shapes seen in Fig. ^ arises mainly from the 
behavior of the resistance, which assumes a linear scal- 
ing with the thickness Nb in the metallic regime and in 
the thick insulating regime (although it has an additional 
constant there, when extrapolated to Nb — Q), but has a 
rapid crossover to the linear regime for the thin insulator 
(semilogarithmic plot shown in the inset to Fig. |l4|). 
One may wish 
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that 



to conclude from Fig. 
lated insulating barriers are superior to metallic barri- 
ers since the parameter ^b/Nb can be reduced to much 
smaller values than in the metallic cases before the char- 
acteristic voltage becomes reduced. But such a view is 
erroneous, because the significantly smaller values of ^f, 
for the insulating barriers means that the barrier thick- 
nesses where the characteristic voltage starts to decrease 
are indeed smaller for the correlated insulator. What 
can be inferred from the figure, however, is that once 
one reaches the critical thickness where the barrier has a 
metal-insulator "transition," the characteristic voltage is 
very strongly dependent on the thickness of the junction. 
Hence variations in the thickness of the barrier can have 
a large effect on the performance of a junction with a 
correlated-insulator barrier. In particular, variations in 
the thickness could make junctions appear to have "pin- 
holes" because slightly thinner areas can have greatly 
enhanced Josephson coupling. This can possibly explain 
why it appears to be more difficult to attain small spreads 
in junction properties for high-Tc-based junctions, even 
if the barrier is pinhole free because the proximity to the 
thickness triggered metal-insulator transition generates 
"intrinsic pinholes" within the correlated insulator. 

It is also interesting to examine how the anomalous 
average behaves as a function of the thickness of the 
barrier as well. We find the following result shown for 
Ufk = 4 in Fig. Il^: once the thickness is larger than 
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FIG. 14. Characteristic voltage plotted versus the inverse 
of the effective thickness of the barrier on a log-log plot. Us- 
ing the correlation length extracted from Fig. hS, allows us to 
plot the characteristic voltage against a measure of the Thou- 
less energy ~ ^-KksT^^^ /N^ . Such a plot should show 
scaling behavior, according to the quasiclassical theory; we 
find this to be approximately true for the metallic junctions 
{Ufk = 2, circles; and Ufk = 4, squares), but the correlated 
insulating barrier has a much sharper dependence on the bar- 
rier thickness (including an "oscillation") and the scaling of 
the quasiclassical theory breaks down. Note the sharp onset 
of insulating behavior at a thickness Nt ~ 7 for Ufk ~ 6. 
Inset is the resistance versus barrier thickness for Ufk ~ 6. 
Note the sharp location of the metal-insulator transition near 
Nt = 7. 

the bulk coherence length (i.e. for all barriers simulated 
with TVf, > 5), we find that the shape of the anomalous 
average is identical for all thicknesses for the planes that 
lie within the superconducting region and that penetrate 
two to three planes into the barrier. What this tells us 
is that the thickness of the barrier is not influencing the 
shape of the anomalous average except within the barrier 
itself, so the oscillations are a property of the bulk su- 
perconductor coming in contact with the barrier. These 
results are also true for the Ufk = 2 and Ufk = 6 cases, 
but we don't show those results here, because the agree- 
ment is essentially the same as seen in the Ufk — 4 figure 
below. 

Finally, we examine the current-phase relation. We 
find similar behavior to that found when one fixed the 
thickness and tuned the correlation strength through the 
metal-insulator transition. For the weakly correlated 
metal {Ufk = 2), we find that the maximum of the 
current-phase relation occurs for phase differences much 
less than tt/2 as expected. As the barrier is made thicker, 
the maximum first overshoots tt/2 and then as Nt, > 20, 
it settles down at 7r/2. This is also seen in the pseu- 
dogap phase Ufk — 4, but the spread in values for the 
maximum of the current-phase relation remains clustered 
closer to 7r/2. By the time we reach the insulator phase 
Ufk = 6, the maximum monotonically increases from 
below 7r/2 for thin barriers to 7r/2 for thick barriers (with 
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FIG. 15. Anomalous average versus plane num- 

ber for {Ufk ~ 4) and thicknesses ranging from 
A'^i,=l,2,5,10,15,20,30, and 40 (the thickness of each barrier 
is obvious from the range of the plots). Note how the shapes 
of these curves are identical for the regions close to the super- 
conductor-barrier interface at planes 25-34 (and within the 
superconducting region to the right or the left for the right 
or left interface, respectively). Since this shape stops chang- 
ing after Nb = 2, we conclude that the oscillations are an 
intrinsic property of the bulk superconductor terminated on 
the barrier. 

no overshoot). 

V. CONCLUSIONS 

In this work we have examined what happens as the 
barrier of a junction is tuned from a metal to an insulator 
for short coherence-length s-wave superconductors. We 
studied the transition both as a function of the correla- 
tion strength and of the barrier thickness. We found a 
number of interesting results. First, in regimes where the 
critical current density approaches that of the bulk super- 
conductor, self-consistency is important in determining 
the current-phase relation, and it is modified dramati- 
cally from simple sinusoidal behavior. As the correlations 
increase, and the current density decreases, the sinusoidal 
behavior is restored, but in some cases, the maximum of 
the current-phase curve overshoots 7r/2 and then becomes 
sinusoidal only at an even larger correlation strength. 
Second, we found that as the barrier becomes more in- 
sulating, the anomalous pair-field average develops oscil- 
lations on the order of the Fermi wavelength, which can 
be quite substantial in amplitude (up to about twice the 
bulk anomalous average). These oscillations are tied to 
the superconductor-insulator interface, and depend little 
on the thickness of the barrier once the thickness is larger 
than about twice the bulk coherence length. Third, we 



14 



found that the critical current has a nontrivial depen- 
dence on the thickness of the barrier — while it decays 
exponentially with thickness, it also has a power-law pref- 
actor that varies with correlation strength, and deviates 
sharply from the quasiclassical prediction. The barrier 
coherence length decreases, of course, as the correlation 
strength increases into the insulating regime. Fourth, we 
found that the characteristic voltage has rich behavior. 
It is maximized for weakly correlated metallic barriers 
for thin junctions and the Ambegaokar-Baratoff result is 
recovered at strong correlations. As the barrier thickness 
increases, the maximum in the metallic region is reduced, 
but the Ambegaokar-Baratoff result fails as the correla- 
tions increase, with the voltage increasing linearly with 
Ufk over a wide range of correlation strengths. The 
intermediate thickness junctions have the most interest- 
ing behavior — the voltage initially decreases, has a sharp 
increase at the metal-insulator transition, and then de- 
creases in the largo correlation limit. Thick insulating 
barriers have very low characteristic voltages and strong 
temperature dependence, as expected at finite tempera- 
tures, since the junction resistance decreases rapidly as 
the temperature is increased in the insulating regime. 
We also saw that self-consistency can renormalize the 
Ambegaokar-Baratoff limit, reducing it by about 10% for 
the single-plane barrier. Fifth, we saw a dramatic devi- 
ation from the quasiclassical predictions as the barrier 
becomes insulating due to strong electron correlations. 
The characteristic voltage remains high for a larger range 
of Thouless energy than in metallic junctions, and then 
decreases very rapidly as the barrier passes through a 
critical thickness where the metal-insulator transition oc- 
curs. This behavior leads to the possibility of "intrinsic 
pinholes." 

This work shows that correlation effects, and the inter- 
play between superconductivity and oscillations brought 
on by the underlying Fermi surface, play increasingly im- 
portant roles in short-coherence- length superconductors. 
In particular, optimizing the characteristic voltage of a 
junction near the metal-insulator transition is possible, 
but requires a careful tuning of the thickness of the bar- 
rier, the proximity to the metal-insulator transition, and 
the operating temperature of the device. Correlated in- 
sulating barriers can mimic effects due to pinholes be- 
cause the Joscphson coupling depends very strongly on 
the thickness leading to an "intrinsic" pinhole effect. In 
the future, we plan on extending this work to d-wave 
superconductors for direct applications to high- Tc super- 
conductors. 



VI. ACKNOWLEDGEMENTS 

We are grateful to the Office of Naval Research for 
funding under grant number N00014-99-1-0328. Real- 
axis analytic continuation calculations were partially sup- 
ported by the National Computational Science Alliance 



under grant number DMR990007N (utilizing the NCSA 
SGI/CRAY ORIGIN 2000) and were partially supported 
by a grant of HPC time from the Arctic Region Super- 
computer Center. We wish to acknowledge useful dis- 
cussions with T. Van Duzer, J. Ketterson, T. Klapwijk, 
J. Luine, J. Mannhart, I. Nevirkovets, N. Newman, J. 
Rowell, and S. Tolpygo. J.K.F. thanks the hospitality of 
the IBM, Almaden Research Center, where this work was 
completed. 



1 B. D. Josephson, Phys. Lett. 1, 251 (1962). 

2 K. K. Likharev, Rev. Mod. Phys. 51, 101 (1979). 

^ For a recent review see the special issue of Superlattices 
and Microstructures 25, No. 5/6 (1999). 
A. Kastalskii, A.W. Kleinsasser, L.H. Greene, R. Bhat, 
P.P. Milliken, and J.P. Harvison, Phys. Rev. Lett. 67, 3026 
(1991); V.T. Petrashov and V.N. Antonov, JETP Lett. 54, 
241 (1991). 

C. W. J. Beenakker, Rev. of Mod. Phys. 69, 731 (1997). 
® C.J. Lambert and R. Raimondi, J. Phys.: Condens. Matter 
10, 901 (1998). 

K.A. Delin and A.W. Kleinsasser, Supercond. Sci. Technol. 

9, 227 (1996); R.S. Decca, H.D. Drew, E. Osquiguil, B. 

Maiorov, and J. Guimpel, Phys. Rev. Lett. 85, 3708 (2000). 
® A. S. Barrera and M. R. Beasley, IEEE Trans. Magnetics 

MAG-23, 866 (1987). 
^ A. Bardas and D. V. Averin, Phys. Rev. B 56, R8518 

(1997). 

A. Levy Yeyati, A. Martm-Rodero, and F.J. Garcia- Vidal, 

Phys. Rev. B 51, 3743 (1995). 

J.C. Cucvas, A. Martm-Rodcro, and A. Levy- Yeyati, Phys. 

Rev. B 54, 7366 (1996). 
12 W. Belzig, F.K. Wilhelm, C. Bruder, G. Schon, and A.D. 

Zaikin, Superlattices and Microstructures 25, 1251 (1999). 
" J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 
1* S.-K. Yip, Superlattices and Microstructures 25, 1213 

(1999). 

" W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 

(1989). 

M. Potthoff and W. Nohing, Phys. Rev. B 59, 2549 (1999). 
" P. Miller and J. K. Freericks, J. Phys.: Conden. Matter (to 
appear). 

1* F. Sols and J. Ferrer, Phys. Rev. B 49, 15913 (1994). 

1^ Li-Fu Chang, S. Chaudhuri, and P.F. Bagwell, Phys. Rev. 

B 54, 9399 (1996). 
2° N. Agrait, J.G. Rodrgio, C. Sirvent, and S. Vieira, Phys. 

Rev. B 46, R5814 (1992). 

21 J. Hubbard, Proc. R. Soc. London, Scr. A 276 238 (1963). 

22 J. Bardeon, L. Cooper, and J. Schrieffer, Phys. Rev. 108, 
1175 (1957). 

2^ P. G. de Gennes, Superconductivity of Metals and Alloys 

(Addison- Wesley, 1966). 
2^ P. Nozioros and S. Schmitt-Rink, J. Low Temp. Phys. 59, 

195 (1985). 



15 



M. Keller, W. Metzner, and U. Schollwock, |cond- 
m |at/010l"047[ 

A.J. Legget in Modern trends in the theory of condensed 
matter, edited by A. Pekalski and J. Przystawa (Springer, 
Berlin, 1980). 

L. M. Falicov and J. C. Kimball, Phys. Rev. Lett. 22, 997 
(1969). 

P. G. J. van Dongen, Phys. Rev. B 45, 2267 (1992); 
P. G. J. van Dongen and C. Leinung, Ann. Phys. (Leipzig) 

6, 45 (1997). 

J. Bardeen, Rev. Mod. Phys. 34, 667 (1962). 
^° Y. Nambu, Phys. Rev. 117, 648 (1960). 
'^^ E. N. Economou, Green's Functions in Quantum Physics, 

App.B (Springer- Verlag, Berlin, 1983). 

A. Khurana, Phys. Rev. Lett. 64 1990 (1990). 

L O. Kulik and A. N. Omelyanchuk, Fiz. Nizk. Temp. 3, 
945 (1977) [Sov. J. Low. Temp. Phys. 3, 459 (1977)]. 
O. N. Dorokhov, JETP Letters 36, 318 (1982). 
K. M. Schep and G. E. W. Bauer, Phys. Rev. Lett. 78, 
3015 (1997). 

Y. Naveh, V. Patel, D.V. Averin, K.K. Likharev, and J.E. 
Lukens, Phys. Rev. Lett. 85, 5404 (2000). 
^"^ V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 486 
(1963). 

O. Gunnarsson and J. E. Han, Nature 405, 1027 (2000). 

B. K. Nikolic and P. B. Allen, Phys. Rev. B 63, R020201 
(2001). 

*° J. Clarke, Proc. R. Soc. London, Ser. A 308, 447 (1969); 
J. G. Shepherd, Proc. R. Soc. London, Ser. A 326, 421 
(1969). 

P. G. de Gennes, Rev. mod. Phys 36, 225 (1964). 

S.-K. Yip, Phys. Rev. B 58, 5803 (1998). 

H. Curtois, Ph. Gandit, D. Mailly, and B. Pannetier, Phys. 

Rev. Lett. 76, 130 (1996). 

H. J. Fink, Phys. Rev. B 14, 1028 (1976). 

A. D. Zaikin and G. F. Zharkov, Sov. J. Low. Temp. Phys. 

7, 184 (1984). 

M. Yu. Kupriyanov and V. Lukichev, Sov. J. Low. Temp. 
Phys. 8, 526 (1982). 
''■^ P. Dubos, H. Courtois, F. K. Wilhelm, A. D. Zaikin, and 
G. Schon, Phys. Rev. B 63, 064502 (2001). 
M. Janssen, Phys. Rep. 295, 1 (1998). 



16 



